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ABSTRACT 
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A group of stations in the North American Arctic region have 
been analyzed for statistical determination of temperatures at 
mandatory -pressure levels. For each station the temperature at a 
key level, peculiar to that station, has been forced in at the 
first step and retained at each subsequent step, in the development 
of stepwise regression equations which predict temperatures at the 
mandatory levels. In general, eight-step predictions, in terms of 
inter-level thicknesses, were found to give optimum specification 
of the desired temperatures. The best estimate of the regional 
atmosphere, which is conditionally dependent upon the existence of 
an extreme 1% probability of the forcing-level temperature, is 
obtained with a high degree of confidence. 
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I . INTRODUCTION 



Although many factors must be considered in determining the 
thermal environment of airborne ordnance, the most dominant are the 
flight conditions and the atmospheric temperature profile, \lhile air- 
craft and missile performance capabilities are well defined, there is 
a great deal of variability in the atmospheric temperature profile 
which may be traversed by the vehicle. Therefore many model atmo- 
spheres have been developed for meteorological purposes, but none is 
completely satisfactory for use in aerothermodynamic heating 
calculations . 

Supplemental standard atmospheres have been developed by 
COESA(1966) to provide information on spatial and temporal variability 
across selected latitude bands and by seasons. However the supplement 
standards are intended to be means with respect to the latitude bands 
and seasons,’ and are not representative of the extremes which must be 
considered for aerodynamic design purposes. Richard and Snelling (1971) 
proposed the adoption of a revised MlL-STD-210 hot and cold model atmo- 
spheres which represent global extreme values at successive 2 km eleva- 
tions. These extremes are deviations from the station mean at the 
elevation involved, which are exceeded at the 1%, 5% and 10% probability 
levels, and were derived through analysis of extreme data from global 
meteorological temperature maps analyzed at constant pressure levels. 
More recently a second draft, MIL-STD-210B, has been published (1972). 
Both of these extreme sets of profiles pose serious design problems when 
changes of elevation are involved. 
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The most important limitation of the two recently proposed 
extreme-atmosphere sets is that they cannot be realistically utilized 
in conjunction with flight profiles and trajectories which involve 
changes in altitude. The temperature extremes at differing altitudes 
are in most cases representative of different geographical locations 
rather than of conditions occurring in real time over a real location. 
Polar and tropical standard and extreme atmospheres have been developed 
for such applications but are still subject to the deficiencies noted 
above. 

TABLE 1. Proposed locations of the 1% world-wide extreme temperatures 
at indicated levels (AFCRL, 1970). 



Level 


Pressure 

altitude 

(ft) 


1% Cold extreme 
location 


1% Warm extreme 
location 


SFC 


0 


Oymyakon, USSR 


Insalah, Algeria 


850 mb 


4780 


bymyakon, USSR 


Insalah, Algeria 


700 mb 


9882 


Hall Beach, NTJT * 


Babylon, Iraq 


500 mb 


18287 


Resolute, Niff - 


New Delhi, India 


300 mb 


30066 


Thule , Greenland* 


New Delhi, India 


200 mb 


38661 


Thule, Greenland* 


Alert, NWT * 


150 mb 


44646 


Karachi , Pakistan 


Alert, NT'JT * 


100 mb 


53084 


Singapore 


Thule , Greenland* 



The object of this study was to develop the most probable verti- 
cal profile, conditional with the existence of a 1% extreme at a given 
level, for each station set forth in Table 1 (but here limited to those 
marked by an asterisk). Each station of Table 1 has been associated 
with a listed 1% extreme for the indicated month according to a pre- 
liminary draft of the MIL-STD-210B proposed atmosphere, which is 
graphically illustrated in Fig. 1. The asterisk-marked stations all 
occur in tVie Canadian Arctic and are either extremely cold at the 
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indicated level in winter or extremely warm in the summer stratosphere. 
This subset of the stations in Table 1 was isolated from the full set 
because the data were readily available in conveniently summarized 
form from checked rawinsonde data. The financial assistance and kind 
cooperation of the Commander, Naval Weather Service is gratefully 
acknowledged in acquiring the data on magnetic tape from the National 
Environmental Data Service, Ashville, North Carolina. For each station 
and month of interest, four years (1967-70) of atmospheric soundings 
were processed to determine the temperature means, standard deviations 
and inter- level correlations required for application of the multiple 
regression analysis procedure to be used in the determination of the 
most probable temperature profile at the 1% probability level. 
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II. DATA PROCESSING 



For each of the stations marked with an asterisk in Table 1, four 
years of checked rawinsonde data, for the period 1967 through 1970, 
was provided by the Naval Weather Service Environmental Detachment, 
Ashville. Those stations listed under the cold extreme of Table 1 
were subject only to January data samples, and those under the warm 
extreme, only to July samples. The individual sounding samples were 
arranged in the sequential format indicated in Table 2 and recorded 
on magnetic tape. Height and temperature data for the indicated 
pressure levels were arranged with six sets in each of the first three 
rows and five sets in the fourth row. Note that the data is arranged 
at 50 mb increments from 1000 to 200 mb inclusively and at somewhat 
smaller increments thereafter. 



TABLE 2. Arrangement of sounding temperature and geopotential 
height data by pressure levels k=l,...,23 in each sounding. 



Index 


k = 


1 


2 


3 


4 


5 


6 


level 


(mb) 


1000 


950 


900 


850 


800 


750 


Index 


k = 


i 


8 


9 


10 


11 


12 


level 


(mb) 


700 


650 


600 


550 


500 


450 


Index 


k = 


13 


14 


15 


16 


17 


18 


level 


(mb) 


400 


350 


300 


250 


200 


175 


Index 


k = 


19 


20 


21 


22 


23 




level 


(mb) 


150 


125 


100 


80 


70 





With the data in the format of Table 2, it is convenient to 
determine the mean temperature T(P ) and standard deviation o for the 

vJ V/ 
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indicated station levels specified by Table 1. The temperature T(P ) 

J 

is to be used in a multiple regression analysis to determine the most 
probable vertical temperature profile at all mandatory-sounding levels 
associated with that ”f orcing’-level temperature" T(P^). A subset of 
such station data is further analyzed for determination of the above 
statistical parameters for those soundings having 10% and 1% extreme 
temperatures at the forcing level. The T^(.IO) and T^(.Ol) extreme 
temperatures can then be forced into the derived full-sample regression 
procedure to obtain the corresponding extreme profiles. 

TABLE 3. Nominal temperature extremes at indicated stations and 
pressure levels at the 10% and 1% values of a normal distribution. 



January Cases 










Full 


Sample 




Forcing 






Sample 


Size 10% 


Station 


Level T^ 


Tj(.lO) 


Tj(.Ol) 


Size 


Extreme 


Hall Beach 


T(700) 


-35.6 


-42.3 


243 


26 


Resolute 


T(500) 


-47.1 


-51.3 


237 


25 


Thule 


T(300) 


-62.0 


-67.3 


213 


13 


Thule 


T(200) 


-63.5 


-70.6 


213 


18 


July Cases 


Thule 


T(IOO) 


-41.3 


-39.8 


236 


23 


Alert 


T(150) 


-39.7 


-37.9 


244 


22 


Alert 


T(200) 


-38.6 


-35.7 


244 


18 


Table 3 


lists the 


forcing level 


pressures 


and sample sizes 


> with 


the 10% and 


1% limiting 


temperatures 


for each station under considera- 


tion. Estimates of the 


temperatures 


at 


the probability 


extremes 
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of 10% and 1% have been computed assuming a Gaussian distribution at 



each J-level. These estimates, Tj(.lO) and Tj(.Ol) of Table 3, are 
derived using the computed full sample-means T(Pj) and standard devia- 
tions Oj of the 1967-70 samples together with the well-know, n 
relationships 



where the plus sign is to be used for a warm extreme and the minus sign 
for a cold extreme. 

Analysis of the data revealed an insufficient number of soundings 
occurring at the 1% extreme of T(Pj) for a test of the statistical 
significance of these results. It was decided therefore, to give 
primary consideration to the 10% extreme T(P^) soundings and to deduce 
the properties of the 1% subset as a limiting case of the 10% sample. 
Thus the same multiple regression treatment used to predict the T(P^) 
mandatory-level temperatures for the full sample has been performed 
upon the "l6% subset" containing T(Pj) extremes. 



T (.10) = T i 1.2817a 

«J «J • 

Tj(.Ol) = T j i 2.3267a 



J 



J 



( 1 ) 
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III. THE REGRESSION METHOD 



A. THE STEPWISE REGRESSION PROCEDURE 

Many problems in research require the extensive analysis of large 
amounts of data, and the process for handling it should be made as 
automatic and rapid as possible. To this end, the Biomedical Computer 
Programs (Dixon, 1966) were developed. One of these programs, BIMED 
02R, provides stepwise regression analysis capability and will be 
utilized here. 

BIMED 02R computes, in a stepwise manner, a sequence of multiple 
linear regression equations. That single variable is added to the 
equation at each step which produces the greatest reduction in the pre- 
vious unexplained variance. This variable is also the one which has 
the highest partial correlation with the dependent variable at the 
particular step in the analysis of the variance. Equivalently it is 
the variable which would give the highest F-statistic upon entry at 
that step. The F-statistic upon entry F^, is expressed at step r as 
(Dixon 1966) 

F (l,n-r-l) - % (C.E.V. ,r) - % (C.E.V. ,r-l) 

^ % (U.E.V.,r) 

where 

% (C.E.V. ,r) is the percent cumulative explained 
variance at step r, 

% (C.E.V.,r-l) is the percent cumulative explained 
variance at step r-1, 

% (U.E.V.,r) is the percent unexplained variance 
remaining at step r. 
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This study utilized a statistical model expressible in the form 



Here in (2) is the temperature predictand at each of the eight 
mandatory levels chosen alternately as 



^ 1 ’ ^ 2 ’ ^ 3 ’ "^ 4 ’ "^ 5 ’ ^ 6 ’ ^ 7 ’ "^ 8 ’ "^9 

but with ^ T(Pj), since the forcing-level temperature is used as the 

forced predictor for each of the eight mandatory level specifications 

for the station under consideration. Thus for each of the stations in 

Table 3, eight of the following nine dependent variables were specified 

for use in the multiple regression procedure: 

T ^5 the temperature at 1000 mb, T(IOOO) 

T^, the temperature at 850 mb, T(850) 

T^, the temperature at 700 mb, T(700) 

T^, the temperature at 500 mb, T(500) 

T^, the temperature at 300 mb, T(300) 

Tg, the temperature at 200 mb, T(200) 

T^, the temperature at 150 mb, T(150) 

Tq, the temperature at 100 mb, T(IOO) 

o 

Tg, the temperature at 70 mb, T(70). 

T(P ) is forced into the specification of each of the eight possible 
j 

predictands in the manner symbolized by Eq. (2), called a struc- 

ture function, represents a linear combination of up to seven inde- 
pendent variables, comprised of layer thickness-values and the derived 
coefficients which meet the requirements of the stepwise least squares 
technique. The thickness variables Zj^, were generated within the BIMED 
02R program by utilization of the so-called transgeneration tecimique 
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which takes the difference of successive tabulated heights in a 
sequential manner from the surface to 70 mb. This procedure produced 
twenty-two independent thickness variables which are directly related 
to the relevant mean virtual temperature in the layer ^k+1^ ’ 

the hypsometric equation. These thickness predictors proved to be 
highly useful in the specification of T^, 

In the specification of the dependent variable T^, a statement is 
permitted as to the maxim\jm number of independent variables to be 
tested for admission to the regression equation. In the tests con- 
sidered here, an upper limit of eight independent variables was found 
in general to give optimal specification of each of the dependent 
variables. Still another convenient feature of the BIMED 02R program 
is that the forced entry into the regression of the temperature value 
being tested as an extreme , the so-called for cing -level temperature, 
is allowed even though its contribution to the explained variance may 
become small at the final step in the regression analysis. 

Tables 4 through 10 list the variables selected and coefficients 
of Eq. (2) for each station extreme level considered. The listed 
structure functions are often highly dependent upon the thickness 
variables nearly adjacent to the level represented by T^^. However, 
diversity of thickness-variable selection does occur. 
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TABLE 4. Full-sample regression equations ^ with associated thickness- 
variables at Hall Beach, NWT, for the specification of the January 
temperatures at indicated levels using 1^(700) as the forcing-level 
temperature. 



h 



^2 

^2 

^4 

^5 

^6 

^6 

^7 

^8 

^8 

^9 

^9 



= -268.4072 + 0.1473T2 + 1.0 

= 0.555Z, + 0.532Z. - 0.521Z, + 0.061Z, - 0.035Z,, + 0.081Z,, 

1 2 3 5 15 16 

-0.059Z^g 

= -285.3345 - 0.0541Tg + 1.0 

= -0.014Z, - 0.016Z_ + 0.319Z- + 0.324Z, - 0.020Z., + 0.008Z- 

1 2 3 4 7 9 

= -293.3438 - 0.0755T_ + 1.0 X, 

3 4 

= 0.106Zg + 0.109 Z^q + O.lOlZ^j^ + 0.095Z^2 ~ 0.036Zj^g 

= -277.0920 + 0.0092Tg + 1.0 X^ 

= -0.041Z,. + 0.171Z,, + 0.094Z_ - 0.014Z., 

13 14 15 16 

= -244.0715 + 0.0905T, + 1.0 X, 

3 6 

= -0.026Z, - 0.017Z., - 0.015Z_ + 0.059Z,, + 0.201Z,^ - 0.026Z,„ 

4 7 15 16 17 18 

= -262.7598 + 0.0306Tg + 1.0 X^ 

= -0.020Z, - 0.037Z,., + 0.l45Zi„ + 0.114Z,_ - 0.018Z„„ 

4 1/ lo iy /U 

= -272.9868 + 0.00651'^ + 1.0 Xg 

= -0.003Zj^ 2 ■ 0.026Zj^g + 0.099Z2Q + 0.105Z2j^ - 0.0U7Z^^ 

= -254.3952 + 0.0447Tg + 1.0 Xg 

= -0.022Z, - 0.008Z„ - 0.005Z_ - 0.029Z„, + 0.308Z.., 

4 9 15 21 22 

is specified in terms of heights z as 



0.555(Zg3Q ^iqqq) + ^-^^^^^'900 ^950^ *^*^^^^^850~^900^ 

+0.061(z^^Q-ZgQg) - 0*035 (z25q-Z3qq) + ^ *081(z2qq-Z25o^ 



0*059(Zj^^q 
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TABLE 5. Full-sample regression equations, with associated thickness- 
variables at Resolute, NWT, for the specification of the January 
temperatures at indicated levels using T^CbQO) as the forcing-level 
temperature. 

T, = -318.9878 - 0.0034T, + 1.0 

1 4 1 

X, = 0.410Z, + 0.426Z. - 0.392Z. + 0.233Z, + O.OIOZ,^ + 0.068Z,^ 

1 1 2 3 4 15 18 

-0.0562^2 

T. = -269.0266 + 0.0200T, + 1.0 X. 

2 4 2 

X. = -0.022Z- - 0.009Z^ + 0.325Z. + 0.284Z, - 0.008Z. 

2 1 2 3 4 6 

T. = -289.0706 = 0.0791T, + 1.0 X. 

3 4 3 

X. = -0.090Z, + 0.169Z. + 0.133Z^ + 0.1242^ + 0.170Z. - O.OllZ^ 

3 4 5 / o 6 9 

= -262.3394 + 0.0421T, + 1.0 X^ 

5 4 5 

X^ = -0.007Z, - 0.011Z_ - 0.0342,. + 0.159Z,, + 0.1002,^ - 0.0172,. 

5 1 12 13 14 15 16 

T. = -272.8584 + 0.0232T, + 1.0 X. 

6 4 6 

X. = -O.OllZ, ^ + 0.064Z,. + 0.2062,^ - 0.0352,. 

D 15 Id 1/ lo 

= -264.4431 + 0.0184T, + 1.0 X_, 

7 4 7 

X^ = -0.0112^2 " 0.0322^^ + 0.133Z^g + 0.1222^^ - 0.0182^^ 

To = -270.4817 + 0.0114T, + 1.0 X. 

8 4 8 

X. = -0.0052^ + 0.0202,. + 0.0832.. + O.IZIZ^, - 0.0582^0 

o / ly ZU Zl 2 . 2 . 

T. = -269.8740 - 0.0072T, + 1.0 X. 

9 4 9 

X. = -0.0052. + 0.0042,, - 0.0052,. - 0.0402^, + 0.3252.. 

9 2 14 16 21 22 
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TABLE 6. Full-sample regression equations, with associated thickness- 
variables at Thule, Greenland, for the specification of the January 
temperatures at indicated levels using TjC300) as the forcing-level 
temperature . 

T^ = -257.2419 - 0.0353T^ 4- 1.0 

X, = 0.534Z, + 0.554Z- - 0.364Z„ - 0.037Z„ - 0.041Zt„ - 0.022Z,, 

1 1 2 3 9 10 16 



+0.043Z 
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T^ = -272.4780 + 0.0045T^ + 1.0 X 2 

X 2 = -O.OI 9 Z 2 + 0.338Z^ + 0.315Z^ - 0.020Z^ - 0.030Z^ 

Tj = -272.8140 + 0.0Q93T^ + 1.0 X^ 

X„ = 0.057Z„ - 0.118Z, + 0.166Z. + 0.162Z, + 0.108Z^ + 0.102Z„ 

3 3 4 5 6 7 8 

T, = -257.5325 + 0.0588T^ + 1.0 X, 

4 5 4 



X, = 0.074Z„ 4- 0.073Z,^ 4- 0.119Z,, 4- 0.125Z,„ - O.Q42Z^„ - 0.019Z 



^4 . 

T, = -246.8040 4- 0.0773T, 4- 1.0 X^ 
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^13 



^14 



X, = -0.008Z, - 0.035Z_ 4- 0.059Z,, 4- 0.198Z. ^ - 0.013Z,_ 

6 1 15 16 17 18 

T^ = -270.2495 4- 0.0083T^ 4- 1.0 X^ 

X^ = -0.009Z^ 4- 0.006Z^ - 0.038Z^^ 4- 0.156Z^g 4- 0.081Z^g 
Tg = -275.9387 4- 0.0118T^ 4- 1.0 Xg 

X- = -0.006Z„ 4- 0.012Z, 4- 0.061Z„_ 4- 0.114Z., - 0.036Z„„ 

8 2 4 20 21 22 

Tg = -283.4604 - 0.0250T^ 4- 1.0 Xg 



X„ = -0.006Z, 4- 0.011Zt„ 4- 0.012Z^, - 0.027Z,. 4- 0.277Z„„ 

9 1 10 16 19 22 
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TABLE 7. Full~sample regression equations, with associated thickness- 
variables at Thule, Greenland, for the specification of the January 
temperatures at indicated levels using T^(200) as the forcing-level 
temperature. 




-331.3064 - 0.2659T, + 1.0 

6 1 

0.526Z^ + 0.5432^ - 0.346Z^ - 0.035Z^ - 0.039Z^^ + 0.086Z^g 
-0. 1567^2 

-272.2056 - 0.0038T, + 1.0 X. 

6 2 

-0.019Z„ + 0.337Z, + 0.314Z, - Q.020Z. - 0.029Z, 

2 3 4 5 6 

-273.0007 + 0.0049T, + 1.0 X„ 

6 3 

0.057Z., - 0.118Z, + 0.165Z. + 0.162Z, + 0.109Z.. + 0.102Z„ 

3 4 5 6 7 8 

-269.6765 + 0.0154T, + 1.0 X, 

6 4 

0.074Z„ + 0.074Z.^ + 0.118Z,. + 0.124Z,„ - 0.045Z,. - 0.0062^, 

9 10 11 12 13 14 

-264.6812 + 0.0508T- + 1.0 X^ 

6 5 

0.009Z^ - 0.055Z^„ + 0.179Z_, + 0.091Z_ - 0.019Z,, 

6 13 14 15 16 

-195.1660 + 0.2897T, + 1.0 X., 

6 7 

O.OIOZ^^ - 0.021Z^g - 0.091Z^^ + 0.161Z^g + O.OSlZ^g 

-283.8765 - 0.02621^ + 1.0 Xg 

0.006Zj^ + 0.008Z^g + O.O 6 OZ 2 Q + 0.114Z2j^ - O.O 36 Z 22 

-280.3738 - 0.0026T, + 1.0 X„ 

6 9 

0.009Zg + O.OlbZj^g - 0.026Z^g + 
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TABLE 8. Full-sample regression equations, with, associated thickness- 
variables at Alert, N\<JT, for the specification of the July tempera- 
tures at indicated levels using TjC200) as the forcing-level 
temperature. 



X, 






X. 



X. 



X. 



-273.6511 - 0.3203T, + 1.0 

6 1 

0.425Z^ + O. 373 Z 2 - O.ZAOZ^ - 0.033Zg - 0.027Zg + 0.029Zj^2 

+0.032Z, , 

16 

-304.2915 - 0.0173Tg + 1.0 X 2 

0.026Z, + 0.024Z„ 4- 0.142Z„ + Q.399Z, - 0.009Z^ - 0.007Z_ 

1 2 3 4 7 15 

+0.036Z^g 

-273.5454 - 0.0034T, + 1.0 X_ 

6 3 

0.034Z^ - 0.132Z, + 0.164Z. + 0.172Z, + 0.127Z^ + 0.114Z_ 

3 4 5 6 7 8 

-0.007Z„ 

-278.9707 +. 0.0340T, + 1.0 X, 

6 4 

0.006Z„ 4- 0.078Z„ 4- 0.082Z,. 4- 0.120Z,, 4- 0.112Z_„ - 0.053Z,„ 

L y lu XX xz xj 

4-O.OlOZ, , 

14 

-248.0954 4- 0.1898T, 4- 1.0 X^ 

6 5 

-O.OllZ^ 4- 0.047Z,, 4- 0,032Z,„ - 0.174Z^„ 4- 0.264Z,, 4- 0.082Z,. 

5 11 13 13 14 15 

-0.040Z, , 

16 

-258.6167 4- 0.0495T, 4- 1.0 X^ 

6 7 

0.015Z, - O.OIOZ^. - 0.004Z.. - 0.047Z, ^ 4- 0.155Z,„ 4- 0.101Z,„ 

4 11 15 1 / lo ly 

- 0.014Z2Q 

-272.3201 4- 0.0112T, 4- 1.0 X„ 

6 8 

-0.006Z- 4- O.OIOZ, - 0.004Z^, - 0.004Z,, + 0.061Z-. 4- 0.134Z„, 

2 4 14 16 20 21 

-O.O 63 Z 22 

-273.5842 4- 0.0138T, 4- 1.0 X„ 

6 9 

-0.006Z^ - O.OIOZ 2 4- 0.022Zg - O.OOSZ^g 4- 0.008Z^g - O.OI 3 Z 22 

4 -O. 27 IZ 22 



23 



TABLE 9. Full-sample regression equation, with associated thickness- 
variables at Alert, Nl'JT, for the specification of the July tempera- 
tures at indicated levels using Tj(150) as the forcing -level 
temperature. 

= -207.5527 - 0.0898T^ + 1.0 

X, = 0.424Z, + 0.374Z^ - 0.248Z. + 0.024Z, - 0.045Z„ - 0.035Z^ 

1 1-2 j D 8 y 

+0.022Z^2 

= -329.3796 - 0.1059T^ + 1.0 X^ 

X^ = 0.026Z^ + 0.0242^ + 0.149Z^ + 0.387Z^ - 0.008Z^^ + 0.055Z^^ 

= -274.1506 - 0.0104T^ + 1.0 X^ 

X„ = 0.033Z^ - 0.127Z, + 0.160Z^ + 0.168Z, + 0.125Z^ + 0.112Z„ 

3 3 4 5 6 7 8 

T, = -271.5181 + 0.0310T^ + 1.0 X, 

4 7 4 

X^ = 0.080Zg + 0.084 Z^q + 0.112Z^^ + 0.1042^^ ~ 0.037Z^^ 

= -276.9827 + 0.1162T^ + 1.0 X^ 

X^ = 0.047Z,, + 0.032Z,„ - 0.175Z._ + 0.264Z, , + 0.0812,^ - 0.0252,, 

5 11 12 13 14 15 16 

T, = -259.7761 + 0.1023T^ + 1.0 X, 

6 7 6 

X, = 0.0072, + 0.0382,, + 0.1842,,, -0.023Z,„ + 0.0162^, 

o 1 Id 1/ 19 21 

To = -285.4023 - 0.0316T,, + 1.0 X„ 

o / o 

Xo = -0.005Z„ + 0.0102. - 0.0022,^ + 0.0622^0 + 0.1362^, - 0.0622^0 

o 2 4 lo 20 21 22 

Tg = -273.5064 + 0.0135T^ + 1.0 X^ 

X- = -0.0072- + 0.0132- - 0.0042,^ + 0.0072, ^ - 0.012Z-, + 0.270Z^o 

y z j JLU ij zi zz 
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TABLE 10. Full-aample regression equation, with associated thickness 
variables at Thule, Greenland, for the specification of the July 
temperatures at indicated levels using TjClOO) as the forcing-level 
temperature . 

T^ = -2i5.4731 + 0.3238Tg +1.0 

Xt = 0.360Z^ + 0.374Z„ - 0.075Z„ - 0.034Z, - 0.041Z, + 0.020Z, ^ 

1 1 -2 J 4 6 15 

-0.021Z, , 

16 

T^ = -303.4197 - 0.0272Tg + 1.0 Xg 

X 2 = 0.047Z^ + O.O 48 Z 2 + 0.105Zg + 0.401Z^ - 0.012Z^ - 0.004Zg 

+0.018Z2^ 

Tg = -319.8106 - 0.1370Tg + 1.0 Xg 

X., = -0.081Z, + 0.164Z^ + 0.156Z, + O.lllZ., + 0.116Z„ + 0.047Z^, 

3 4 5 6 7 8 21 

-O.O 3 IZ 22 

T, = -276.6890 + 0.0313T„ + 1.0 X, 

4 8 4 

X, = O.OllZ., + 0.092Z„ + 0.094Z,_ + 0.096Z,, + 0.096Z,„ - 0.033Z,., 

4 2 9 10 11 12 13 

T^ = -266.3704 + 0.0764T„ + 1.0 X^ 

D O _) 

X^ = 0.014Z., - 0.073Z,- + 0.214Z,. + 0.070Z,_ - 0.016Z-,, 

5 7 13 14 15 16 

T^ = -255.2715 + 0.1766T„ + 1.0 X, 

6 8 6 

X, = 0.019Z^, - 0.035Z_ + 0.061Z^, + 0.123Z,., + 0.103Z,o - 0.044Z„„ 

6 14 15 16 17 18 20 

T., = -184.3870 + 0.0909T_ + 1.0 X, 

X., = -0.030Z, - 0.050Z„ + 0.031Z,, - 0.014Z_ + 0.068Z,., + 0.008Z,„ 

7 6 9 11 16 17 18 

+0.02lZ^g 

To = -257.1123 - 0.020T- + 1.0 X„ 
y o y 

Xo = O.OllZ, + 0.002 Z,t - 0.008Z^„ - 0.003Z,. + 0.24lZ„o 

9 4 11 13 14 22 
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B. RELATED STATISTICAL PARAMETERS 



Available as outputs from the BIMED 02R program are several related 
statistical parameters such as: 

a, multiple R^ after r steps 

b, standard error of estimate, S.E. 

c, mean value, T 

d, standard deviation, a 

e, F-statistic after r steps, 

2 

f, percent reduced variance, R^. 

2 

The significance of R^ in relation to S.E. is given by (Crow et al., 
1955) as 

(S.E.)^ = a^((N-l)/(N-r-D) (1-R^) 

where 

2 ^ - 2 

a = Z (T.-T) / N. (3) 

i=l ^ 

It should be noted that S.E. is a measure of the error in specification 
of T^^ after application of the regression equation estimator of 
form (2). 

2 

In (3) , o is the variance of the temperature sample and 

N is the sample size 

i is the sample- element index 

r is the number of predictors selected 

R is the multiple correlation coefficient after r steps. It 

2 

follows that the percent unexplained variance (1 - R ) , for sample sizes 
of 213 to 244, as contained here, is closely approximated by 

(S.E. / a)^ = 1 - R^. 
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The percent explained variance R , also called the percent reduced 
variance, is then given by 

= 1 - (S.E. / a)^. (4) 

2 

The mean value of R , found in the prediction of the 1000 mb tempera- 
tures using the full-data set, was 0.8478 and 0.7480 with the 10% 

2 

extreme-data set. For all levels above 1000 mb, the pooled mean R for 
the full-data set was 0.9672 and 0.9596 for the 10% extreme-data set. 

The scores at 1000 mb are somewhat lower than those of the upper levels 
and this is to be expected, in that the surface variabilities bring 
periodic and random-effect terms into the specification of the 
temperature at 1000 mb. 
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TABLE 11, Regression statistics at mandatory-pressure levels at Hall Beach, NWT, using T(700) as the 
forcing level. Part (a) refers to the full-data January sample; part (b) refers to the nominal 10% 
extreme sample of T(700) ; part (c) lists the sounding level mean temperatures corresponding to 
the 1% extreme cold set of T(700) . 
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TABLE 12. Regression statistics at mandatory levels at Resolute, NWT, using T(500) as the forcing-level 
temperature. Part (a) refers to the full-data January sample; part (b) refers to the nominal 10% cold 
extreme sample of T(500) ; part (c) lists the sounding-level mean temperatures corresponding to the 1% 
extreme cold set of T(500) . 
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TABLE 13. Regression statistics at mandatory-pressure levels at Thule, Greenland, using T(300) as the 
forcing-level January temperature. Part (a) refers to the full-data sample; part (b) refers to the 
nominal 10% cold extreme sample of T(300); part (c) lists the sounding-level mean temperatures corre- 
sponding to the 1% extreme cold set of T(300) . 




30 



TABLE 14. Regression statistics at mandatory-pressure levels at Thule, Greenland, using T(200) as the 
forcing-level January temperature. Part (a) refers to the full-data sample; part (b) refers to the 
nominal 10% cold extreme sample of T(200); part (c) lists the sounding- level mean temperatures corre- 
sponding to the 1% extreme cold set of T(200). 
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TABLE 15. Regression statistics at mandatory-pressure levels at Alert, NWT, using T(200) as the forcing- 
level July temperature. Part (a) refers to the full-data sample; part (b) refers to the nominal 10% cold 
extreme sample of T(200) ; part (c) lists the sounding-level mean temperatures corresponding to the 1% 
extreme cold set of T(200) . 
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TABLE 16. Regression statistics at mandatory-pressure levels at Alert, NWT, using T(,150) as the forcing- 
level July temperature. Part (a) refers to the full-data sample; part (b) refers to the nominal 10% cold 
extreme sample of T(150); part (c) lists the sounding-level mean temperatures corresponding to the 1% 
extreme cold set of T(150). 
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TABLE 17. Regression statistics at mandatory-pressure levels at Thule, Greenland, using TClOO) as the 
forcing-level July temperature. Part (a) refers to the full-data sample; part (b) refers to the nomina 
10% cold extreme sample of T(IOO); part (c) lists the sounding— level mean temperatures corresponding to 
the 1% extreme cold set of T(IOO). 
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IV. ANALYTICAL PROCEDURE USING THE REGRESSION TFXHNIQUE 



A. RESULTS OF THE REGRESSIONS 

Part (a) of Tables 11 through 17 lists the mean, standard deviation, 
multiple correlation coefficient and standard error of estimate for the 
full-data sample of each station and at all mandatory levels except the 
forcing level. These statistics were computed after application of the 
appropriate regression equations listed in Tables 4 through 10. Part 
(b) of Tables 11 through 17 gives similar results, after application of 
the identical regression equations used in part (a) , to a regression 
analysis of the data set which contains only the ”10% extreme” set of 
forcing -level sample atmospheres. This means that f(P^) was chosen if 
it lay within the particular 10% Gaussian distribution extreme under 
consideration. The remainder of the mandatory levels of the extreme 
sounding, T^(.IO), were then subjected to the same regression analysis 
as the full-data samples of part (a). This nominal set of extreme 
atmospheres v?as selected by choosing those soundings for which 
T(P^) ^ T^(.IO) for the warm extremes, and T(P^) _< T^(.IO) for the cold 
extremes. These extreme 10% values T^(.IO) are defined by the Gaussian 
distribution limits listed in Eq. (1). 

Part (c) lists the mean temperature profile defined at mandatory 
levels which correspond to the set of soundings having T(P ) at the 1% 

J 

extreme. There were two ways of selecting the 1% extreme cases: (i) 

application of an a priori test to determine the values T(P ) which 

J 

fell within the required range indicated by the Gaussian distribution; 
(ii) if there were no values of l(Pj) fitting this constraint, then 
either the two or three coldest (for cold extremes) or warmest values, 



35 



were chosen a posteriori from the regression-determined 1% set as repre- 
senting the set of soundings corresponding to the 1% extreme T(Pj)- 
values. The corresponding set of T(P^) values were then simply averaged 
and considered to represent the expected vertical profile having T(P^) 
in the 1% extreme. 



B. THE TEST FOR SIGNIFICANT SOUNDING DIFFERENCES 

In view of the small sample sizes in parts (b) and (c) of 
Tables 11,..., 17, it was not feasible to test for significant differ- 
ence by random sampling of the 10% set relative to the 1% set of T(Pj) 
generated soundings. It was decided therefore to test for significant 
differences in defining the T^ values predicted in each part (b) of 
Tables 11,..., 17 against those T^ values of part (a) obtained b)^ the 
identical prediction process, but using the full-data sample: that is, 

in case (a) there were no _a priori limiting conditions set up regarding 
a particular value of T(Pj). 

In order to check the latter hypothesis, a Student t-test was 
performed to test whether the 10% samples of prediction means of part (b) 
V7ere significantly different from those of the full-data sample means 
of part (a) . The pooled t-statistic relating to the T^^ prediction 
differences of the part (a) and part (b) means may be formulated after 
Brooks and Carruthers (1953) with (N+n-(r+2)) degrees of freedom, after 
computing the two sets of predictive mean values and their standard 
errors. The t-statistic is formulated as follows; 



t(N+n-(r+2)) = 



^"^FULL 



(N-9)a^ + (n-2)o 



E, .10 



(N+n-r-2) 



i. + 1. 

N n 






(5) 
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Here the T’s are the part (a) and the part (b) sample temperature means 

(at the mandatory levels) based upon identical prediction equations* 

a is the part (a) standard error of estimate, and o m the part 
E t , • J.U 

(b) standard error of estimate. The number of degrees of freedom asso- 
ciated with the two values are (N-r-1) and (n-2), respectively, in 
parts (a) and (b) with 

N = the full-sample size 
n = the nominal 10% extreme-sample size 
r = the number of predictors, r = 8. 

The pooled t-statistic may be computed from the horizontally adjacent 
sections of Tables 11,..., 17 and with the proper number of (N+n-(r+2)) 
degrees of freedom for t. Tables 18 and 19 of Appendix 1 list the 
resulting values at mandatory levels for each station under consideration. 

These t-values were tested against critical t ^^-values (e.g. , 

Crow et al, 1960) with the same number of degrees of freedom used in 
(5) . This number of degrees of freedom usually fell into the tabular 
class for t-values corresponding to an infinite number of degrees of 
freedom. Thus for example t ^ , with an infinite number of degrees of 
freedom, is t = 2.326. Tables 18 and 19 show one value belov^z the 99% 
confidence limit specified by t = 2.326. That value, corresponding 
to the 1000 mb level of Thule (with = 300 mb) , still yields a con- 
fidence limit of 92%. The random chance of the T^ specifications of 

parts (a) and (b) of Tables 11,..., 17 being drawn from the same popu- 

12 

lation if t = 7.0 is approximately one chance in 5 X 10 and is much 

smaller for t > 7.0. Most of the t-values in Tables 18 and 19 

exceeded t = 7.0 with exceptions principally at P = 1000 mb. 
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C. THE TYPE-(C) PROFILE CLASSIFICATION 

As noted previously, it was not feasible to employ a t-test utiliz- 
ing the extremely small sample size of part (c) of Tables 11,..., 17. 
However, it is useful to explain physically the pattern of trends in 
the sample means T^ of part (c) relative to those of (b) . The code 
symbol ”D” for ’’direct”, used in part (c) , means that is dis- 
placed directly in the same sense from as is from 

the level under consideration. The code symbol ”M” for ’’mixed” means 
that 1 t^(c) - T^(a) I < I that a consistent displace- 

ment in the temperature extreme does not occur in progressing across the 

columns of T -values. This code symbolism was used for all cases con- 
M 

sidered and is summarized under part (c) of Tables 11,..., 17. To 

some extent, the nature of these displacements, whether dix^ect or mixed, 

is indicated also by Fig, 2,..., Fig. 8, which show only the comparison 

of the T^. (c) and T. . (a) profiles. 

M M 



38 



I 



V. RESULTS SUMMARY FOR INDIVIDUAL STATIONS 



A. WINTER CASES 

This section will consider those cases investigated as cold ex- 
tremes and summarized in Tables 11,..., 14 and Figs. 2,..., 5. 

1. Hall Beach, M^^r; T(Pj) = Tj(700) 

Table 11(a) shows the results of the stepwise regression 
applied to the specification of the January full-sample temperatures 

(T^ ” *^1» *^2^ *^4^ *^5^ ^6^ *^7^ *^8 mandatory levels of 

Hall Beach by an equation of form (2) . An eight-predictor equation for 
Tj^, at each level, was generated by use of the BIMED 02R program with 
the results shov7n. Very high multiple correlation coefficients were 
obtained, indicating a functional relationship for T^^ at all levels. 
Table 11(b) lists analagous results for the data sample corresponding 
to the 10% cold extreme at the forcing level, with T^j(b) being deter- 
mined by the same specification equation used in part (a) . Part (c) 
presents the resultant estimated mandatory- level means derived 

from the five cases which fell into the 1% extreme forcing level sub- 
set of part (b) . Part (c) then, gives the expected vertical profile 
corresponding to the 1% extreme at *^(Pj) = T(700). It is suggested 
that the mean *^^j(c) profile, relative to that of T^(b), should be phys- 
ically as well as statistically explainable from the conditions imposed 
by the choices of P^.. At Hall Beach, the latter selection at 700 mb 
was consistent with the existence of an extreme-cold troposphere, 

T^(c) < T^^(b), and an extreme-warm stratosphere for winter. The only 
case of a mixed classification M, occurs for T^j(500), v/hcre according 
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FIG. 2. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Hall Beach, NOT, corresponding 
to the set of 1% cold extreme occurrences of Tj(700). The thin 
solid line depicts the January mean Hall Beach vertical sounding 
(1967-70) . 



AO 



Pressure (m 



to Fig. 2, the crossover between the mean and extreme soundings occurs. 
The nature of the profile is then consistent with the displacement 

of the profiles except at the crossover point near the 

trop opause. 

2. Resolute, NInTT; T(P^) = T^(500) 

In the case of Resolute, Table 12(c) shows a somewhat different 
result from that discussed above for Hall Beach. At Resolute x^ 7 e now 
start with the premise of an initially specified very cold 
T (c) < T (500). This level is excluded from being an effective 

J • U JL 

crossover point in the results of Fig. 3 except in the region of P=300 mb 
where a slight tendency tox^ard a warm-stratosphere does occur. Compar- 
ing the two cases, the selection of an extreme cold T^(500) dictates the 
presence of a cold stratosphere coupled with a cold troposphere below. 

For nearly all levels at Resolute, the T^,(c) - T,,(a) trend bears out 
that indicated by T^,j(t>) “ sense of maintaining the T^^(c) 

atmosphere relatively cold at virtually all levels. In contrast v^7ith 
the situation at Hall Beach, the Resolute type (c) atmosphere more 
nearly resembles the case of a cold tropospheric vortex with an over- 
lying stratospheric cold vortex, Xi?hereas the Hall Beach type (c) case 
seems to occur with the warm stratospheric anticyclone aloft. 

3. Thule, Greenland; T(P ) = T (300) 

j a 

This case is graphically depicted in Fig. 4, which shows that 
the A priori condition of a 1% cold temperature extreme at P^=300 mb at 
Thule is associated with a cold stratosphere as well as a cold tropo- 
sphere at all levels above and below Pj=300 mb. The extreme vertical 
profile for this regime (T(P^ = T(300)) at Thule resembles that at 
Resolute, where there was also evidence of a similar coupling betv^ocn 
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FIG. 3. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Resolute, N\\'T, corresponding to 
the set of 1% cold extreme occurrences of T^(500). The thin solid 
line depicts the January mean Resolute sounding (1967-70) . 
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FIG. 4. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the 1% cold extreme occurrences of 1^(300). The thin solid line 
depicts the January mean Thule sounding (1967-70). 
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FIG. 5. The heavy solid line shovzs the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the 1% cold extreme occurrences of T (200) . The thin solid line 
depicts the January mean Thule sounding" (1967-70) . 
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cold stratospheric temperatures in case (c) and a cold-extreme upper 
troposphere. Now, however y Fig. 4 shows no crossover point between the 
mean and extreme profiles. Likewise Table 13 shows the direct transi- 
tion D-type classification applicable at all levels of the mean 
temperature profile . 

4. Thule, Greenland; T(P^) = 1^(200) 

This case differs from the previous case (T(P ) = T (300)), in 

j j 

that the latter was associated by regression development with T^^(200) = 
-62.83^C, whereas the coldest 1% input temperature T^(200) was -74.70^C, 
The imposition of the latter condition into the regression analysis at 
Thule (Tj(200)) is associated with the coldest stratospheric regime 
encountered during the entire study of this paper. With T^(200) = 
-74.70°C, all regression-generated stratospheric temperatures are typi- 
cal of those in the cold Arctic vortex at and above 200 mb. However 
Fig. 5, shows that in this case there is a crossover in the T^^(a) , T^j(c) 
profiles. Furthermore by comparison of Table 14(c) with Tables 11(c), 
12(c) and 13(c), it is clear that this case corresponds to the warmest 
of the four winter tropospheres discussed here. 

5. Combined Winter Regimes, Summary 

The analysis of the four sets of type (c) winter extremes has 
indicated that the 300 mb level is a key indicator in the type of atmo- 
sphere to be expected in the Canadian Arctic region. In summary, 
these regimes have been identified by the four stratifications: 



(a) 


Hall Beach 


Tj(700), 


cold 


extreme 


(b) 


Resolute 


Tj(500), 


cold 


extreme 


(c) 


Thule 


Tj(300) , 


cold 


extreme 


(d) 


Thule 


Tj(200), 


cold 


ex Lr erne 



45 



Case (a) corresponds to an abnormally cold troposphere with a 
resultant abnormally warm stratosphere (Fig. 2). Cases (b) and (c) 
correspond jointly to cold stratospheres and cold tropospheres (Figs. 3 
and 4), but neither cold layer exists at a particular extreme. Finally 
case (d) is characterized by a record cold stratosphere and a near- 
record warm lower troposphere. Cases (a) and (d) show negative corre- 
lations of temperature across the tropopause, a property documented at 
mid-latitude stations by Cole and Nee (1965) . However cases (b) and 
(c) have a consistent positive correlation throughout the atmosphere 
and have yet to be documented fully. 

B . SUMMER CASES 

This section will consider those cases as having warm extremes at 
forcing levels and are summarized in Tables 15, 16 and 17 and Figs. 6, 

7 and 8. The same stepwise regression procedure was applied to the 
July data of the warm extreme stations as was applied to the January 
data of the cold extreme stations. The results are presented below. 

1. Alert, N\^fT; T(P^) = T^(200) 

Table 15 lists the statistical results of the analysis of Alert 
with T(P ) = T (200. Consistently high multiple correlations are 

J J 

shown indicating a high degree of confidence in the predicted T^^ values 
including those of part (c) . This confidence is reinforced by the con- 
sistent direct transition D-type classification applicable at most 
levels of the mean temperature profile Tj^j(c) . Mixed tendencies occur 
at 100 mb and 70 mb but are characterized by extremely small tempera- 
ture differentials (less than 0.5^C) in the trends. A 

crossover of the extreme and profiles exists between 500 mb 




FIG. 6. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Alert, NWT, corresponding to the 
1% warm extreme occurrence of T (200) . The thin solid line depicts 
the July mean Alert sounding (1967-70) . 
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pressure (m 



and 300 mb (Fig. 6), indicating an extreme warm stratosphere overlying 
a relatively cool troposphere. This phenomenon is also characteristic 
of mid-latitude negative inter-level temperature correlations sho™ by 
Cole and Nee (1965) to exist near the level of the tropopause during 
periods of cyclone-anticyclone variability. 

2. Alert, NWT; T(P^) = 1^(150) 

Table 16 corresponds to the analysis of Alert with T(Pj) = 

T (150) . Very high multiple correlation coefficients in both parts (a) 

J 

and (b) are again shown yielding a resultant high degree of confidence 
in the diagnostic procedure of isolating "extreme" atmospheres. 

Figure 7 indicates that a marked similarity exists between the T^^(a) 
and T (c) profile pairs and the corresponding pair for Alert of T (200) . 

1 V I 

Table 16 indicates that the T.,(c) regression-generated atmosphere is 

M 

displaced primarily in a "direct" sense relative to the T^^(b) atmo- 
sphere. It should also be noted that the forcing extreme T^(150) at 
the nominal 10% and 1% values differ only by about l^C . The comparison 
between generally of the same character as that of 

Tj^^(c) relative to T^^(a) , (Fig. 7). 

The same conditions of warm stratosphere and relatively cool 
troposphere, as well as the crossover between 500 mb and 300 mb, sug- 
gest a similar comparison with regard to this climatic regime as with 
the preceding one (Alert Tj(200)). 

3. Thule, Greenland; T(Fj) = T^(IOO) 

The final table of regression statistics, Table 17, corresponds 
to the analysis of Thule at T(P^) = Tj(lOO). Once again there are very 
high multiple correlation coefficients in both parts (a) and (b) . In 
addition the atmospheric profiles, both in the mean and 1% extreme 
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Temp (^C ) 



FIG. 7. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Alert, NWT, corresponding to the 
1% warm extreme occurrence of T (150). The thin solid line depicts 
the July mean Alert sounding (1967-70) . 



A9 




FIG, 8. The heavy solid line shows the mean regression-determined 
vertical temperature sounding over Thule, Greenland, corresponding 
to the 1% warm extreme occurrence of T (100) . The thin solid line 
depicts the July mean Thule sounding (1967-70) . 
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(Fig, 8), were nearly identical to the corresponding profile pairs in 
the previous two warm cases (Figs. 6 and 7). This indicates a strong 
tendency of the synoptic processes acting at Alert in July to also be 
acting at Thule during the same period. 

4. Combined Summer Regimes, Summary 

Each of the July cases is based upon the assumption of a wam 
extreme forcing- level temperature existing at or above 200 mb. In each 
case, there proved to be a concurrent warm temperature at 300 mb, that 
is, stratospheric warming above 200 mb has extended downward to 300 mb 
in July. At tropospheric levels, however, the mean temperatures of 
parts. (b) and (c) are generally colder than the full-sample means. 
Figures 6, 7 and 8 are graphical presentations of the mean temperatures 
shown in parts (a) and (c) of Tables 15, 16 and 17. 

Three factors should be recognized as being relevant to the 
July data considered. (i) In two of the three cases, the nominal 1% 
extreme temperature at the forcing level T(P^) did not occur, and the 
class was redefined by a. posteriori methods to comprise the set of two 
warmest observations at the forcing level. (ii) The difference between 
the parts (b) and (c) mean mandatory level profiles (Tj^(b) and T^j(c)), 
at and above 300 mb, never exceeded l.O^C. (iii) The dominant feature 
of all of the summer case-studies was the consistent occurrence of 
stratospheric- tropospheric temperature reversals noted in the 
preceding paragraph of this section. 
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VI. CONCLUSIONS 



The multiple regression procedure presented for development of 
vertically consistent model atmospheres is statistically dependent upon 
the input condition of a known extreme temperature at a given flight 
level. The procedure has been shown to be highly definitive in deter- 
mining the most probable vertical temperature profile over the stations 
of interest. Extremely warm stratospheric temperatures are not asso- 
ciated with extremely warm tropospheric temperatures, thus integrated 
ballistic density anomalies tend to be nearly compensated. Likewise, 
extremely cold tropospheric temperatures tend to be reversed in the 
stratosphere. These ti'jo points become extremely important in the 
determination of ballistic parameters. As a result, aerothermodynamic 
calculations can be made for flight conditions involving changes in 
altitude by using an atmospheric temperature profile which is consis- 
tent with the proposed MIL-STD-210B value at the primary level of 
interest, but is more realistic and less severe at all other levels. 
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TABLE 18. Student's t-scores for non-forcing mandatory levels at indicated North American Arctic cold 
extreme stations as computed from the full data samples and 10% extreme data samples using equation (5). 
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TABLE 19. Student's t-scores for non-forcing mandatory levels at indicated North American Arctic warm 
extreme stations as computed from the full data samples and 10% extreme data samples using equation (5).. 
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